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Filing Date 


22 May 2000 
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PRELIMINARY AMENDMENT AND 
REQUEST NOT TO ENTER ANY AMENDMENTS 
UNDER ARTICLE 34 OF THE PCT 



Box PCT 

U.S. Patent and Trademark Office 
P.O. Box 2327 
Arlington, VA 22202 

Sir: 

Prior to calculation of the filing fee for this 
application, please enter the following Preliminary 
Amendment : 

If any amendments have been filed under Article 34 
of the Patent Cooperation Treaty, applicant requests that 
these amendments not be entered as a part of this national 
filing. 

IN THE ABSTRACT 
Please add an abstract as appearing on the 
separate page herewith. 
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IN THE SPECIFICATION 

Page 1, after the title and between lines 1 and 2 insert 
the following captions and paragraph: 

—CROSS REFERENCE TO RELATED APPLICATION 
— This application is Section 371 application based on 
International Application No. PCT/FROO/01402 filed May 22, 
2000, and published as WO 00/74387A1 on December 7, 2000 
not in English. 

— FIELD OF THE INVENTION — 

Page 1, between lines 5 and 6, insert the following 
caption: 

—BACKGROUND OF THE INVENTION-- 

Page 4, between lines 3 and 4, insert the following 
caption" 

— SUMMARY OF THE INVENTION — 

Page 7, between lines 25 and 26, insert the following 
caption: 

—BRIEF DESCRIPTION OF THE DRAWINGS — 

Page 8, between lines 25 and 26^ insert the following 
caption: 

— DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS— 



IN THE CLAIMS 

Please amend claims 3-4, 6-14, 17, and 19 as follows: 



3. (Amended) Encoding method according to claim 1, 
characterized in that said mesh units are triangles (81) . 



4. (Amended) Encoding method according to claim 2, 
characterized in that said specific encoding comprises the 
following steps: 

- the association, with each of said triangles (81) to 
be encoded, of a square matrix (84) representing said 
triangles to be encoded, by means of a first reversible 
transformation (72) ; 

- the application of a second reversible decorrelation 
transformation (74) to each of said square matrices (84) 
delivering transformed matrices. 

6. (Amended) Encoding method according to claim 25, 
characterized in that said second transformation (74) 
belongs to the group comprising: 



7. (Amended) Encoding method according to claim 16, 
characterized in that said specific encoding comprises a 
step of quantification (75) and of encoding the data of the 
lower part of said transformed matrix. 

8. (Amended) Encoding method according to claim 17, 
wherein at least at least one of the following parameters 
can be modified: 



- the Karhunen Loeve transformation (KLT) , 

- the discrete Fourier transformation (DFT) , 
" the discrete cosine transformation (DCT), 

- and the Walsh Hadamard transformation (WHT) . 



- said first threshold; 



- said second threshold; 

- the structure of said nested meshing; 




- a scaling factor a enabling an expansion or 
compression of said square matrix in the space domain; 

- the type of qualification; 

- the quantification pitch, 

9. (Amended) Encoding method according to claim 1, 
characterized in that said piece of error information 
corresponds to a difference in luminance, representing a 
squared error or an absolute error between a source 
triangle and the corresponding approximate angle. 

10. (Amended) Encoding method according to claim 1, 
characterized in that, for each mesh unit, the following 
choice is made: 

- if said piece of error information is below said 
first threshold, the hierarchical division for the mesh 
unit considered is interrupted (132); 

- if said piece of error information is above said 
first threshold, but below said second threshold, the 
application of said hierarchical division (134) is 
continued; 

- if said difference in luminance is above said second 
threshold, said specific encoding (135) is implemented. 

11. (Amended) Encoding method according to claim 110, 
characterized in that, for all meshing levels other than 
the first level, the data taken into account for the 
encoding are quantified differential values, each 
corresponding to the quantified difference between an 
optimized value for the hierarchical level considered of 



the image to be encoded and a value interpolated from the 
meshing of the lower hierarchical level. 

12. (Amended) Encoding method according to claim 1, 
comprising a step for optimizing the position of the peaks 
of each of said meshings, minimizing a predetermined error 
criterion. 

13. (Amended) Encoding method according to claim 1, 
comprising a step for defining a quaternary subdivision 
tree used to identify the mesh units using said specific 
encoding and the mesh units whose contents are interpolated 
from the values at the nodes defining said mesh units. 

14. (Amended) Encoding method according to claim 1, 
comprising a step for matching said mesh units with the 
contents of the image, implementing inversions of diagonals 
of the quadrilaterals formed by all the pairs of adjacent 
triangles. 

17. (Amended) Decoding method according to claim 15, 
characterized in that said specific encoding comprises the 
following steps: 

- the association, with each of said source triangles 
to be processed, of a square matrix representing said 
source triangle, by means of a first reversible 
transformation; 

- the application of a second reversible decorrelation 
transformation to each of said square matrices delivering 
transformed matrices. 



characterized in that said complementary decoding 
advantageously comprises the following steps of 
reconstruction : 

a) the application of an inverse transformation to 
said second reversible transformation on said transformed 
matrices, delivering said reconstructed square matrices; 

b) the association, with each of said reconstructed 
square matrices, of a corresponding reconstructed triangle 
by means of an inverse affine transformation of said first 
reversible transformation; 

c) the reconstruction of said minimum partition from 
said reconstructed triangle. 

19. (Amended) Decoding method according to claim 15, 
characterized in that said preliminary decoding provides 
for the reading, in the bit stream^ of at least one of the 
pieces of information belonging to the group comprising: 

the number of levels of the hierarchy; 

the identification of the encoding technique used for 
each of the triangles; 

the succession of the differential values of the 
components associated with the nodes of the hierarchical 
meshing; 

the identification of the arcs on which a diagonal 
inversion is made. 



REMARKS 

The above amendments are made to eliminate multiple 
dependent claims, and to cancel omnibus claims. 

Also, please do not enter the amended pages 
submitted under Article 34. This Preliminary Amendment 
conforms the claims to those considered in the International 
Preliminary Examination Report. 

Entry of the Amendment prior to calculating the 
filing fee is respectfully requested. 

Respectfully submitted. 
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Robert M. Angus, Reg. No. 24,383 
WESTMAN, CHT^MPLIN & KELLY, P. A. 
Suite 1600 - International Centre 
900 Second Avenue South 
Minneapolis, Minnesota 55402-3319 
Telephone : 612/334-3222 
Telefacsimile: 612/334-3312 



MARKED -UP VERSION OF REPLACEMENT PARAGRAPHS 



Page 1, after the title and between lines 1 and 2 insert 
the following captions and paragraph: 

--CROSS REFERENCE TO RELATED APPLICATION 
— This application is Section 371 application based on 
International Application No. PCT/FROO/01402 filed May 22, 
2000, and published as WO 00/74387A1 on December 7, 2000 
not in English. 

FIELD OF THE INVENTION— 

Page 1, between lines 5 and 6, insert the following 
caption: 

--BACKGROUND OF THE INVENTION— 

Page 4, between lines 3 and 4, insert the following 
caption" 

— SUMMARY OF THE INVENTION — 

Page 7, between lines 25 and 25, insert the following 
caption : 

—BRIEF DESCRIPTION OF THE DRAWINGS— 



Page 8, between lines 25 and 26, insert the following 
caption: 

--DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS-- 



MARKED-UP VERSION OF REPLACEMENT CLAIMS 



Please amend claims 3-4, 6-14, 11, and 19 as follows: 

3. (Amended) Encoding method according to an y of t ke- 
claim-& 1 and 2 , characterized in that said mesh units are 
triangles (81) . 

4. (Amended) Encoding method according to e-Bry — 

claime 2 — a^^-e a, characterized in that said specific 

encoding comprises the following steps: 

- the association, with each of said triangles (81) to 
be encoded, of a square matrix (84) representing said 
triangles to be encoded, by means of a first reversible 
transformation (72) ; 

- the application of a second reversible decorrelation 
transformation (74) to each of said square matrices (84) 
delivering transformed matrices. 

6. (Amended) Encoding method according to e-^y of tho 

claime 2 tro 5, characterized in that said second 

transformation (74) belongs to the group comprising: 

- the Karhunen Loeve transformation (KLT) , 

- the discrete Fourier transformation (DFT), 

- the discrete cosine transformation (DOT), 

- and the Walsh Hadamard transformation (WHT) . 

7. (Amended) Encoding method according to amy of t-h e 

claims 1 to 6, characterized in that said specific encoding 
comprises a step of quantification (75) and of encoding the 
data of the lower part of said transformed matrix. 



8. (Amended) Encoding method according to any of th e 
claim-s- 1 — t^e — 7, wherein at least at least one of the 
following parameters can be modified: 

- said first threshold; 

- said second threshold; 

- the structure of said nested meshing; 

- a scaling factor a enabling an expansion or 
compression of said square matrix in the space domain™j_ 

- the type of qualification; 

- the quantification pitch. 

9. (Amended) Encoding method according to^ any of the 
claimiS- l-™te — 8-, characterized in that said piece of error 
information corresponds to a difference in luminance, 
representing a squared error or an absolute error between a 
source triangle and the corresponding approximate angle, 

10. (Amended) Encoding method according to irfty of th e 
claime l—fe-e — 9^, characterized in that, for each mesh unit, 
the following choice is made: 

- if said piece of error information is below said 
first threshold, the hierarchical division for the mesh 
unit considered is interrupted (132); 

- if said piece of error information is above said 
first threshold, but below said second threshold, the 
application of said hierarchical division (134) is 
continued; 

- if said difference in luminance is above said second 
threshold, said specific encoding (135) is implemented. 



11. (Amended) Encoding method according to a ny of fefee 
claim-s- 1 — fee — 10, characterized in that, for all meshing 
levels other than the first level, the data taken into 
account for the encoding are quantified differential 
values, each corresponding to the quantified difference 
between an optimized value for the hierarchical level 
considered of the image to be encoded and a value 
interpolated from the meshing of the lower hierarchical 
level . 

12. (Amended) Encoding method according to any of the 
claims 1 — ^ — ^r, comprising a step for optimizing the 
position of the peaks of each of said meshings, minimizing 
a predetermined error criterion. 

13. (Amended) Encoding method according to eey of trie 
claims 1 t o 12, comprising a step for defining a quaternary 
subdivision tree used to identify the mesh units using said 
specific encoding and the mesh units whose contents are 
interpolated from the values at the nodes defining said 
mesh units. 

14. (Amended) Encoding method according to any of the 
claims 1 to 13, comprising a step for matching said mesh 
units with the contents of the image, implementing 
inversions of diagonals of the quadrilaterals formed by all 
the pairs of adjacent triangles. 

17. (Amended) Decoding method according to — e# 

the claime 15 and lo, characterized in that said specific 

encoding comprises the following steps: 



- the association, with each of said source triangles 
to be processed, of a square matrix representing said 
source triangle, by means of a first reversible 
transformation; 

- the application of a second reversible decorrelation 
transformation to each of said square matrices delivering 
transformed matrices, 

characterized in that said complementary decoding 
advantageously comprises the following steps of 
reconstruction : 

a) the application of ean inverse transformation to 
said second reversible transformation on said transformed 
matrices, delivering said reconstructed square matrices; 

b) the association, with each of said reconstructed 
square matrices, of a corresponding reconstructed triangle 
by means of an inverse affine transformation of said first 
reversible transformation; 

c) the reconstruction of said minimum partition from 
said reconstructed triangle. 

19. (Amended) Decoding method according to any of the 
claims 15— -B-e — 18-, characterized in that said preliminary 
decoding provides for the reading, in the bit stream, of at 
least one of the pieces of information belonging to the 
group comprising: 

the number of levels of the hierarchy; 

the identification of the encoding technique used for 
each of the triangles; 

the succession of the differential values of the 
components associated with the nodes of the hierarchical 
meshing^ 
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the identification of the arcs on which a diagonal 
inversion is made. 



09/980108 




VARIABI^E TKANSFOKMATIQIL HIERARCHTCAT. TMAGE ENCODING 

The field of the invention is that of the encoding of still or moving images. 
More specifically, the invention relates to image compression techniques or image 
sequence compression techniques based on the implementation of invertible 
mathematical transformations. 

There are very many kno-wn techniques of image compression to reduce tlic 
quantity of data needed to represent an image or a sequence of moving images,. 
Thus, it is sought especially to reduce the hit rates of the digital signals in order to 
transmit tiiem and/or to store them on a data carrier. 

The invention can be applied especially but not exclusively to the 
transmission of image signals at low bit rates as well as to transmisaion without bit 
rate guarantee, as in the case of transmission made according to the IP ("Internet 
Protocol"). 

Among the many known image-encoding methods, it is possible to 
distinguish especially the ISO- JPEG and ISO-MPEG techniques which have given 
rise to a standard. These encoding techniques rely especially on the implementation 
of transformations, enabling the efficient elimination of redundancy in an image. 

The niEiiri traiisrorxnalioni> used in imagtj compression are; 
\ - the Karhunen Loeve transformation (KLT), 
the discrete Fourier transformation (DFT), 
the discrete cosine transformation (DCT), 
and the Walsh-Hadamard transformation (WHT). 
To oYcroomc these various drawbacks, another method has boon proposed 
m the patent document FR-98 12525 filed on behalf of the same Applicant as in the 
present patent application, entitled "Precede de codage d'images fixes ou animees 
avec reduction et adaptation du debit" "Method for the encoding of still or moving 
images with reduction and matching of the bit rate". 

The object of this technique is a method for encoding a digital image 
designed to produce a bit stream representing this image, the length of the bit stream 
being a function of The desired representation. This method comprises the following; 
steps: 

the defining, in a field of the image to be encoded, of a hierarchical 
meshing comprising a plurality of nested meshings for which the 
peaks of the mesh imits may be pixels of said image; 



r 
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the achieving of optimization of luminance, chrominancej and 
positions on each meshing level; 

the determining, for each mesh unit of said hierarchicai meshing, of a 
difference in luiwinajoce between the image to be encoded and an 
5 interpolated image obtained from the peaks of the nested meshing to 

which the mesh unit considered belongs; and 

the introducing, into the bit stream, of the values (advantageously 
encoded differentially with respect to the previous hierarchical level) 
of positions, luminance and chrotninance of the pealcs of the mesh 
10 units whose luminance difference is greater than a threshold 

difference. 

It will be noted that this technique is not limited to the luminance and. 

chrominance signals but can be appUed to any model of colors. According to this 

technique, and as illustmted in Figure 1, at the end of the meshing step, a quaternary 
15 tree structure 1 1 is built This structm-e is associated with the hierarchical meshing 

12 to manipulate the values (colors and positions) of the peal<s of the mesh units. 

The tree has a number of nodes equal to the number of triangles in the corresponding 

meshing level. Each node 13 of the tree pertains to a single triangle 14 of the 

hierarchical meshing 12, 
20 Once this tree 1 1 is built, a selection is made of the data of the tree to be 

introduced into the bit stream representing the image that will be transmitted and/or 

stored. This selection depends on the desired quantity. 

To mjikc iliiis tjclctjtiunp fur cudti triuriglc a. diffcxciiuc in luixuiifcuiL;c between 

the image to be encoded and the interpolated image is calculated on the basis of th^ 
25 peaks of the nested meshing to which the considered mesh unit belongs. This 

difference is then compared with a threshold difference for each triangle. The value 

of the threshold difference is a function of the desired quality of representation. 

Then, the part of the tree relating to the triangles whose difference in 

luminance is greater than the threshold difference is introduced into the bit stream. 
30 This process of selecting the nodes of the tree by taking an in-depth route is 

illustrated in Figure 2, Only the nodes that are above the borderline 21 are preserved.. 

The threshold difference therefore enables the transmission of the data 

relating to the image as a function of the local quality of these different triangular 

partitions. Indeed, on a textured part, the transmission of the data takes place up to. 
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the last meshing level (the finest meshing) ^d, for the smooUier parts, a rough level* 
proves to be sufficient. 

This technique is also described by Wang et al in the docmnent "USE OF 
TWO-DIMENSIONAL DEFORMABLE MESH STRUCTURES FOR VIDEO 
5 CODING (XP000641037). This document proposes a method for the generation of a 
quadrangular liicrai-chical meshing. An error criterion is used to decide whether or 

not it is necessary to continue the hierarchical division for a given mesh unit 

To obtain adequate quality, when a portion of the image is highly textured, 
it is necessary to plai) for a large number of levels. In this case, the efficiency of the 
10 hierarchical encoding is low. In other words, the hierarchical encoding is well suited 
to relatively simple images but not to the images having highly textured parts. 

In practice, it is not rare for at least a part of the image to be highly textured. 
I his aspect truly limits the pmctical value of the hierarchical encoding as described 
here above. 

15 In particular, it is a goal of the invention to overcome these different 

drawbacks of the prior art. 

More specifically, it is a goal of the invention to provide an image encoding 
method that makes use of the advantages of the technique herein known as 

hierarchical encodings without showing any of its drawbacks. 
20 In particular^, it is a goal of the invention to provide art encoding method of 

this kind that is efficient even when the image to be processed has highly textured 
portions. 

Another goal of the invention, naturally, is that the modifications made in 
the hierarchical encoding technique as described here above should not modify its 
25 efficiency when the image or image portion to be processed has low texturing. 

It is also a goal of the invention to provide a method of this kind that is 
simple to implement and especially does not need any manual intervention (for 
example for the selection of the te^itured zones). 

More generally^ a goal of the invention is to provide an image encoding 
30 method with a good ratio between the quantity of data to be transmitted and/or stored 
and quality of image reconstruction. 

Another goal of the invention is to provide a corresponding decoding 
method that enables the reconstruction of images simply and at low cost (in terms of 
processing time, storage capacity, etc.). 
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It is also a goal of the invention to provide a decoding method of this kind 
with high qualities of scalability (i.e, allowing a rough image to be swiftly 
reconstructed and then gradually refined up to the desired level of refinement). 

These goals ?»s w«11 ;3s nfhers that, sha.11 appear mnre clearly hereinafter are 
5 achieved according to the invention by means of an image encoding method 
comprising the following i>tepi>: 

- the definition of a hierarchical meshing having at least two nested meshing 
levels formed by mesh units defined by peaks (which may be pixels of said image to 
be encoded) 

10 - the determining, for each of said mesh units, of a piece of error 

information between said image to be encoded and ein interpolated image obtEiincd 
from the peaks of the mesh units belonging to the meshing level of the mesh units 
considered; 

- pvTtting a stop to the refining of the mesh units having a piece of error 
15 information below a predetermined first threshold; 

- the implementation of a specific encoding for the mesh units having a 
piece of error information above a second predetermined threshold; 

- the continuation of the refining of the nxesh units having a piece of error 
information above said first predetermined threshold and below said second 

20 predetermined threshold. 

Thus according to the invention, two distinct encoding modes are used 
selectively. The hierarchical encoding or nested meshing encoding is the main 
ciicudiag or the basic encoding, but it is not used systematically: for the image 
portions that require it, another more efficient encoding is tised, 

25 The decision to use this specific encoding is taken hy comparison with a 

threshold. More specificallyj three possibilities are offered at each node considered 
of the meshing: stopping the refining or the division of the mesh units (the quality 
attained for the corresponding image portion being sufficient), going to the next 
hierarchical level in preserving the same hierarchical encoding or xising another 

30 encoding (in particular an encoding better suited to the highly textured portions). 

According to a preferred embodiment of the invention, said specific 
encoding implements at least one reversible transformation (for example a DCT); 
Indeed^ this technique proves to be well suited to the encoding of the tet^tiireA 

portions^ 



5 ; 

Advantageously, said mesh units are triangles (obtained by a partitioning j 
that 13 regular or adapted to the image contents). It will be noted that the use of ai 
transformation such as a DCT on triangles is novel per se . ; 

Thus, said specific encoding may advantageously comprise the follov^nngj 
5 slieps: I 

- the association, with each of said tiianglcs to be encoded^ of a squaiel 
matrix representing said triangles to be encoded, by means of a first reversible' 
transfonnation; I 

i 

- the application of a second reversible decorrelation transformation to eachj 
10 of said square matrices delivering transformed matrices. 

Preferably^ said step of associating a square matrix comprises the following' 
Steps: i 

- the aftine transformation of a source triangle into an isosceles rectangular 
triangle called a reference triangle; 

1 5 - the creation of a square matrix whose lower part includes data representing: 

sftid isosce1.e.<j rftCt?5ngn1fiT triangle^ | 

- the symmetrizing of said square matrix i 
Said second transformation may belong in particular to the group comprising: ' 

- the Karhunen Lo^ve transformation (KLT), i 
20 - the discrete Fourier transformation (DFT), 

- the discrete cosine transformation (DCT), i 

- and the Walsh-Hadamard transformation (WHT). 

According to a preferred embodiment, at least one of the following: 
parameters can be modified: 
25 - said first threshold; i 

- said second threshold; \ 

- the structure of said nested meshing; 

- a scale factor a enabling an expansion or compression of s?*irl square: j 
matrix in the ^pace domain, i 

3t> - the type of quantification; \ 

- the quantification pitch. ! 
Advan.tageoiisly, said piece of error information corresponds to a 

difference in luminance, representing a squared error or an absolute error 
between a source triangle and the corresponding approximate triangle, \ 
35 For each mesh unitj, advantageously the following choice is made: i 

i 

I 
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- if said piece of error information is below said first threshold, the i 
hierarchical division for the mesh unit considered is interrupted; \ 

- if said piece of error information is above said first threshold, but below said; 
second threshold, the application of said hierarchical divi^iion is continued; | 

- if said difference in luminance is above said second threshold, said specific | 
encoding is implemented. | 

The adaptation to a signal component other than tine chrominance is direct. j 
Preferably, for all the meshing levels other than the first level, tlie data taken | 
into account for the encoding are quantified differential values^ each corresponding | 
to the quantified difference between an optimized value for the hierarchical level ! 
considered of the image to be encoded and a value interpolated from the meshing of | 
the lower hierarchical leveL i 
This greatly reduces the quantity of the intbrmation to be encoded. | 
Advantageously, the method of the invention comprises a preliminary step fotj 
optimizing the position of the peaks of each of said meshings, minimizing a | 
predetermined error criterion, | 
In other words, the partitioning is adapted to the contents of the image. j 
Preferably, the method comprises a step for defining a quaternary subdivision 
tree used to identify the mesh units using said specific encoding and the mesh units 
whose contents are inteirpolated from values at the nodes defining said m,esh imits. 

Advantageously, it also planned to have a step for matching said mesh units • 
with the contents of the iraage^ implementing inversions of diagonals of the | 
quadrilaterals formed by all the pausi uf adjacent triangles. \ 
The invention also relates to a method for the decoding of data encoded j 
according to the encoding method described here above. A decoding method of this j 
kind comprises especially the following steps: | 

- the preUminary decoding of said data encoded according to said encoding, | 
implementing a hierarchical meshing^ enabling the deecription of an initial j 
representation; i 

- the complementary decoding of said data encoded by means of said specificj 
encoding, enabling the refining of said initial representation. j 

Thus, the image can be reconstructed gradually, as and when reception and ! 
decoding takes place, in an increasingly refined way (first of all from successive i 
hierarchical levels and then from complementary data of the specific encodmg). It is 
therefore possible to match the quality of the reconstrucLioxi uf tlie hnagc to die ! 



available means (transmission bit rates and processing capaciLy uf aie terniiiml), and | 
even to define several levels of decoding quality (it is possible for a simplified i 
version to have no complementaiy decoding). i 
When Bald specific encoding jmplements a DCT applied to triangles, asj 
described here above, said complementary decoding advantageously comprises fhej 
following steps of reconstruction: i 

- the application of a inverse transformation to said second reversiblej 
transformation on said transformed matrices, delivering saidi 
reconstructed square matrices; I 

- the association, with each of said reconstructed square matrices, of a| 
curresponding reconstructed triangle by means of an inverse afftn^ 
transformation of said first reversible transfortnation; ! 

- the reconstruction of said minimum partition from said reconstructed 
triangle. \ 

Advantageously, said square matrices are recreated from data of a received 
bit stream, whose decoded data are the coefficients of the triangle to b^ 
reconstructed, which form the lower part of said matrix. i 

Preferably, said prelimin,ary decoding provides for the reading, in the biit 
Stream, of at least one of the pieces of information belonging to the group 
comprising: : 

the number of levels of the hi erarchy ; | 

the identification of tibe encoding technique used for each of the triangles; | 

the succession of the differential values of tlic components associated wilh thp 

I 

nodes of said hierarchical meshing ! 
the identification of tlie arcs on which a diagonal inversion is made. | 
Other features and advantages of the invention shaU appear more clearly from 
the following description of a preferred embodiment given by way of a simpje 
illustrative and non-restiictive example, and the appended drawings, of which: | 
Figures 1 and 2 illustrate the general principle of the encoding method 
that is the object of a patent application FR-9S 12525 discussed in tiic 
introduction. More precisely: ! 
Figure 1 shows the correspondence between a nested mesh ai^d 
a quaternary tree; \ 
Figure 2 exemplifies selections of nodes of the tree |of 

i 

Figure 1: j 



! 

i 

Figure 3 is a flow oKart showing the general principle of The method j 

of the invention; j 

Figure 4 referred to in Appendix 1 illustrates the determining of the | 

differential values to be introduced into the bit stream after | 

quantification^ which will be transmitted or stored when thej 

hierarchical encoding is implemented; j 

Figure 5, which is commented on in Appendix 2, illustrates an| 

operation for optimizing the positions of the peaks of a meshing level! 

mid it5i propagation at the higher meshing level; i 

Figures 6A and 6B commented upon in Appendix 3 respectively showj 

a case where the diagonal inversion is excluded and a case where it isj 

I 

possible; j 
Figure 7 is a simplified flow chart of a mode of implementing thej 
specific encoding according to the invention, relying oi^' 
the implementation of a DCT; j 
Figure 8 illustrates the principle of the second and third steps of thq 
method of Figxire 7; i 
Figure 9 is a more precise extract of Figure 9; | 
Figi-tres 10 an.d 1 1 show two modes of quantification that can be used 
in the method of Figure 7 ; i 
Figure 12 illustrates the ?;igzag course of the encoding s^ep of the 

method of Figure 7; ; 
Figure 13 is a simplified flow chart illusirating the choice of the 
processing to be made when the method of the invention and ^ 
hierarchical encoding are implemented in an associated way. j 

1. Basic principles of the i nvention I 
The invention therefore relates to an hnprovement of the technique of nested 

ro.eshings (or hierarchical exicoding techniques) as described in the intrnduction. 

According to the present invention, it is planned indeed, to use a second type of 

encoding for the parts of the image for which this technique of nested nic^liin^s 

proves to be inefficient. | 

As explamed here above, a definition is made first of all, according to ih}s 
nested uiesliing techrdque, of a hierarchical meshing in the field of the image to be 
encoded comprising a pluraUty of nested tneshings. The peaks of these meshes m^y 
be pixels of the image to be encoded. This meshing is for example obtained by 
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regular and successive divisions of the meshes of the rough meshing-j 
Advantageously, the position of the nodes defining the mesh units is adapted to thej 

contents of the image. ; 

The nested meshes technique takes position at a level n (hetween the firsli 
and last mesh level) of meshing* the interpolated image is computed by thc^ 
hierarchical meshing technique and an error image corresponding to the difference iri 
luminance between the original image and the interpolated image is deduce4 
therefrom, I 

Then, the tree pertaining to the n first mesh levels is constructed and the 
luminance difference is computed for each of the triangles of the mesh of the error 
image and a threshold difference S is chosen. The criterion of the difference in 
luminance on a triangle T corresponds to the following mean quadratic error: ! 

! 



Where I is the error image between the mterpolated image and the original 
image on the triangle | 

According to the present invention, for each node, it is determined whether 
the approximation procedure must be stopped, whether it is necessarj^ to continue tl|e 
subdivision of the meshing by refined interpolation with the hierarchical meshiijg 
technique or whether it is necessary to use another encoding technique which p 
advantageously a technique based on reversible Lranbrormation as described here 
below. 

Figure 3 illustrates the general operation of the invention. It starts ^vith the 
25 hierarchical encoding 31 at the first level of the hierarchy. A particular example pf 
implementation is described in detail here below (§ 2). ! 

For each triangle cunsiaexed, a decision 32 is taken for the rest of t^e 

procedure: | 
interruption 33 of the encoding; j 
30 - passage 34 to the following hierarchical level; 

implementation of another type of encoding 3 5 . 
Paragraph 3 shows a mode of implementation of an encoding 35 usin^ a 
DCT transform applied to triangles. Paragraph 4 gives a detailed view of 
exemplary decision mechanism 32. 
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2- Reminders on the nested mesh techniques i 

Here below, we shall specify certain advantageous aspects of the 
implementation of the step 32. The general principles of this technique have alreadyj 
been described in the introduction with rf^ference to Figures 1 and 2. | 

2.1 Variable threshold diiference I 
It is advantageously planned Ihut tlie tliresliold difference will dimlnishj 

during the method to provide for scalability of the encoding. I 
Thus, at each hierarchical level, the thresholding step is reiterated several 

times with increasingly small threshold differences so thatj at each iteration, a group^' 

I 

of additional bits is generated to improve the quaUty of the depiction. j 
To build a final bit stream without redundancy, provision is made fot the 
use of a reference table in which each node of the tree is associated with a value 0 oj; 
I indicating whether the values of the node in question, have aheady been introduced 
into the bit stream. The values of the table are updated as the values of the nodes ax^ 
introduced into the bit stream. Thus, when the threshold difference dhninishes, it ik 
ascertained that the node values to be introduced into the bit stream have not ^^trftady 
been introduced. This enables the introduction into the bit stream of only that part of 
the additional tree that corresponds to the lowering of the threshold difference. ! 

It can thus be decided to use only the first part of the bit stream obtained 
with the highest threshold difference when there is no need for high representation 
quality. To obtain higher representation quality, it is necessary to use a j^eater par^ 
of the bit stream. | 

2.2 Quantification and compression \ 
It can also advantageously be planned to quantify and compress the valuep 

of the peaks before introducing them into the bit stream in order to limit the size 6t 
the stream. The quantification is chosen to be preferably non-uniform. It may fqr 
example rely on the statistical distribution of the luminance values (and chrominance 
and position valucti rc£>pectivcly), j 

At the end of the quantification step, the quantified differential values ai^e 
subjected to a data compression operation designed to limit the length of the bit 
stream* Advantageously, the compression is done by ati adaptive arithmetic encoder 
which has high compression rates in the present case. i 

The statistical distributions of the luminEince, chrominance and positionihg 
values being different from one another, it is preferable to encode these values 
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separately. Thus, the adaptive arithmetic encoder i2> designed to process these; 
different types of values separately, | 

2.3 Differential encoding \ 

i 

In order to reduce the sjTie of the bit stream, it is planned according to one; 
5 improvement to introduce, into the bit stream, not the precise luminancej 
chrominance and posiLiun values of the peaks of mesh units of the nested mes^ihirig 
but differential values, except for the basic meshing. Each differential value 
represents the difference between the optimized value of the hierarchical levejl 
considered of the image to be encoded and an interpolated value obtained from exact 
10 values corresponding to the neighboring peaks of the nested tneshing of the lowei: 

liierai'diical IcvcL I 

An exemplary determining of differential values is discussed in Appendix 1[ 

2.4 Meshing adapted to fM conrenis of the ima^/e j 
The structure of the meshing chosen here above is got by regular 

15 construction. It can therefore sometimes appear to be unsuitable for depicting an 
image that has dJsparities in its contents axid that mixes, for example, uniform regions 
with more textured zones requiring a higher density of peaks. An optimizing of the 
position of the peaks of the lucsliing will shift th.o meshing peak concentrations 
towards the zones that requhe it. i 

20 According to a preferred embodiment^ it is planned also lo adapt Lib 

meshing (position and values Y, U, V of the peaks) to the contents of the image so as 
to exploit its encoding potential to titte best possible extent. This adapting can be 
done at three levels: j 
the optimizing of the position of the peaks of the meshing: the 

25 position of the nodes is modified tts a function of their efficiency in 

representing the image locally; 

the optimizing of the luminance and chrominance values of the peaks 
of the meshing: the values Y, V are optimized to represent the 
original image to the best possible extent; I 

30 - the optimizing of the topology of the meshingt the topology of ijhe 

meshing is modified by means of an operation of diagonal inversion 
so as to further improve the capacity of the meshing to represent the 
image locally. l 
The most immediate visual effect of such an optimization is that the peaks 

35 of tlic meshing come closer to the physical contours of the objects of the image. ! A 
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technique of this kind is presented for example in tlie patent application FR-98 12525 
filed on behalf of the Applicant of the present patent application. 

This operation is done level by level, beginning vs^ith the levelj 
corresponding to the basic Tnei?hmg (level 0). The result of the optimization at thd 
5 level n is then transmitted to the level n+1 and the position of the addition^ peaks of 
the level n+1 is oblaiued from the new position of the level n peaks. This optimizine 
of the position at each level and its propagation to the different levels of the tree ar^ 

illustrated in Figure 5 and discussed in Appendix 2. | 

i 

2.5 pia^onal inversion j 
10 Accordmg to a final improvement, it is also possible to plan for improvmg 

the typology of the meshing. Indeed, it happens that the local structure of th^ 
meshing is not suited to particular features of the image. Although it would be 
possible to refine the meshing until an adapted hicraicliical meshing is obtained, it 
sometimes proves to be simpler to cany out a diagonal inversion operation. | 
1 5 This diagonal inversion operation is done on the last mesh level and consis^ 

in inverting the diagonal in the convex quadrilaterals formed by two adjacent 
triangles of said meshing level if the meshing thus modified provides better quality 
of image recousuniction. | 
This operation of diagonal reversal is presented in Appendix 3 wit^ 
20 reference to Figures 6 A and 6B. I 
2.6J Prnhlp.m ryf highly structured zones j 
Despite these different improvements, it is observed that this technique p 

not sufEcicntly effective in processing certain image portions, especially when the^ 
are highly textured. ] 
25 3. TiCrr encodiTig adapted tn a trlan K"^ar partition | 
The invention therefore proposes the implementation of a transformation, 
for example a DCT tiansformation, adapted to a triangular partition for the image 
portions that require it. The choice of tliese portions is discussed in paragraph 4. 
3,1 Principle 

30 Figure 7 is a general flow chart iUustrating the correspondin£ method. 

The processing with node according to the invention is therefore t^e 
following: 

defhiition 71, on the field of the image to be encoded, of a triangu 
partition which may be adapted to the contents on the field of *e 
35 image portion to be encoded; 
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13 ; 

the determiniag for each element of the paiUtion obtained of th.e\ 
transformations used to associate a reference triangle (72) and then a; 
square (namely a matrix) with each triangular embodiment (73); i 

i 

the performance of a DCT 74 on each of these matilces; ; 
5 - the application of a qnantificalion method 75 and encoding method 7^ 

that may be identical with those of current staxidards, | 
The partitioning 71 may of course be adapted to the contents of the image 
for example by implementing the operations described in paragraph 2.4. I 

10 3,2.1 Affine transformation in an isosceles triang le | 

The second mid third stops 72 and 73 of tha method of the invention are 
illustrated in Figure 8. i 
For each triangular element 8 1 of the partitJonj the upeiation determines thje 
^^ffine transformation 82 used for the association with each unspecified triangle 81 of 
15 a reference triangle 83 which is an isosceles triangle. The reference triangle is then 
cunvcrtcd into a square and more precisely a square matrix S4 by symmetrization 8^, 
More precisely, the first transformation 82 is used to determine the afFme 
transformation by which it is possible Lo pass from any triangle &1 to the reference 
triangle 83 as shown in Figure 9, { 
20 The invertible affine transformation F such that = F(Qj) with Pj = (xj, ^ j) 

£ind Qi = (Xi, Yj) is written as follows: | 

This aftine transformatioii is invciiiblc because the determinant of tjhe 
25 matrix is equal (apart from the sign) to 2A (where A represents the area of the 
unspecified triangle 81), which is assumed to be non-zero. This inverse aTfinc 
transformation is therefore written as follows: | 

(•^3--^?>(3^2-)^l)+(^2-^lX>'3-y3) ! 

30 I 
3,2.2 DCT \ 
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if l(x,yj is the luminance of the pixel having coOTclinaLcb (x,y) and if it is 
assumed that the image to be encoded has been partitioned into M x N sized blocks, 
the application of a block-oriented transfonnation a(x, y, m, n) will produce an image 

F wifh: 

M-IN-L ; 

Fim,n)^^^I(x,y)a(x,y,m,n) (1) 

j;=0y=O 

for ICx,y) = I(x,y) owing to the synimetrization 85. 

Starting with the transformatioa a(x, y, ra, n), a reverse transformation b(x, 
10 y, m, n) can he defined in order to reconstruct the original image 1: 

!(x^y) = X X nm,n)bix,y,m,n) (2) 



The DCT is obtained by assiming: 

15 

, 2c(m)c{n) X2^+l)g7z , ^^ + 1)?^ . 




Witii: 

if w - 0 
else 

20 The second transformation 73, S6 consists of the transposition of the 

information contained in each triangle with an area A in the lower part of a square 
matrix G of iS'(axV2x/l) lines where E represents the higher integer part of the 
value between parentheses and a e R"^^* represents a scale factor that acts on the 
visual depiction of the image, in obtainh:ig an expansion (a>l) or a compression 
25 (Oc<l) in the space domain. 

According to the formulae (1) and (3) we have: 
F(m,n) = F(n, m) 

After the symmetrizatinn of its transformation 74 according to the 
equation (1) generates a matrix that is also symmetrical H. 
30 Consequently^ the information contained in llic lower pait of each matrix G 

being identical to the upper part 75, the use of the block-based DCT transformation 
74 may be implemented, for example by MPEG or JPEG techniques. 
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After the transtbrmation 74, only the lower parts of the matrices H will be 
quantified 75 and encoded 76. 

3.3 Optimization of the encodinO' 

In order to optiini^ft the performance characteristics in terms of the cost of. 
encoding 76, two action means can be implemented:, modulated for example as a 
jfunction of the relevance of die te;ctm-c niadcriying the triangles considered. These" 
are: 

the scale factor a (then a<l is taken); 

the choice of the quantification and especially the amplitude of the: 
quantification pitch chosen. 
3r3.I Quantiftcatio n 

Among the quantifications 75 possible, we can use particularly: 
a uniform quantification; 
a quantification with a zigzag route; 

a quantification by use of a weighting matrix pre-assessed on a 
psycbn-visnal criterion. 
The quantification with a zigzag route consists in initializing the 
quantification process at a value Q^aC= Nvhicli, during the route, at each trace-back, is 
incremented by a value Aac, as illustrated by the arrow 101 of Figure 1 0. 

An exemplary weighting matrix pre-assessed on a psycho-visual criterion 
tba -TPEG standard QM matrix illustrated in Figure 11 or MPEG4. The matrices G 
and QM xnaybe of different sizes. An interpolation will be made of the matrix QM 
bringing thiy uiatrix to the size of G as for JPEG and it is then possible to define a 
quality factor qf acting as a multiplier of the matrix QM, 
3.3J Encoding 

Tlie effective encoding 76 is achieved for example by performing an RLE 
(Run Length Encode) and entropic type of encodmg on a zigzag route 121 shown in 
Figure 12. \ 
4. Combination of the two encoding techniques 

Ihe implemeniation accurding to the invention of the two encoding 
techniques described here above may rely on the method shown in Figui'e 13 . If, f^r 
the given level n, the difference in luminance on a triangle V of the mesh is: 

131: below a threshold difference: the pait of the image inteipolated 
on this triangle has good visual quality and the procedure is stopped 
(132); 
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133 : higher than the threshold difference but below k x S, with k > 1 : 
the method of continuous approximation with the hierarchical mesh 
teclmiquc (134), the part of the interpolated image corresponding to 
an im^igft with average texturing; 

135; greater than k x S with k k 1; the triangle is processed by a DCT 
applied to llic tiianglc of the error image (136)^ ; 
This selection is justified as follows. It is known that: j 



whence: 



therefore: 



\F(m,ni<^^\Kx,y)a(r.y.m,n)\ d aprts CD 

\Fim,n)\^2xY^\lix,y)\ 



It is Hierefore seen that the coefGcient F(m, n) tends towards zero When the 
luminance difference tends towards zero. A low quadratic error leads to coefficients 
AC after low amplitude traasform having lugh chances of being cancelled after 
quantification. • 

Thus, it proves to be more judicious to carry out an affme interpolation, 
which is less csostly than a DCT transformation, on such mesh units. : 

Tbe overall method therefore consists in processing a part of (be hnage by 
the hierarchical meshmg technique and m processing Ibe highly textured parts of Ibis 
image by a DCT according to the present invention applied to triangles of the 

corresponding error image. i 
Here therefore, on the textured part of the error image, a DCT is applied jtn 

the triangles for which the luminance difference is great. ■ 

It must be noted that the transfono.ation-based technique is only an exami|le. 
The invention can use many other complementary encoding techniques, such as for 

example: ; 

methods based on fractal decomposition; the principle of image 
compression in gray levels by the IFS method, also called fi:a9tal 
compression, relies on the expression of tbe contents of Hic image by 
means of the contents themselves. ; 
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It can be seen as a self-quantification of the image. The formalizing; 
of this method arises especially out of the work of Hutchiuson in 1981; 
md the work of Bradley^ Demoko and otb,er research workers at the; 
Georgia Institute of Technology between 1985 and 1988. The first 
automatic algoritlim applying these ideas to image compression was 
proposed by Jacquin in 19S9. j 
Improvements in this technique are proposed in the patent application 
FR-99 00656 entitled "Procede et dispositif de codage k base de 
schemas TFS, a fonctions de collage oscillanteSj proc6de de codage; 
fonction de collage, support de doiin^es et applications 
conrespoadaiits" (Encoding method and device based on fF5? schemes* 
with oscillating collage functions^ encoding method, collage function, 
data medium and corresponding applications)- ; 
methods known as "matching pursuit" methods described for example 
in the article by Ralph Neff and Avideh Zakhor, "Very Low Bit Rate 
Video Coding based on Matching Pursuits", published in IEEE 
Transactions on Circuits and Systems for Video Technology- 
The encoding (of the residue) by matching pm-suit is an iterative 
method that uses a dictionary of redundant functions. At each 
iteration, a search is made for the function that best represents the 
residue obtained in the previous step. Thus, the image is broken down 
on a sequence of atoms that represents it optimally; 
ttie SADCT (Shiipe Adaptive DCT) described for eXJimple by T. 
Sikora and B, Makai in "Shape Adaptive DCT for Generic Coding'^ 
(IEEE Transactions on Circuits and Systems for Video Tet;luiology, 5 
(1 ), pp. 59-62, February 1 995). 
It can also be planned that several of these techniques will be available, the 

selection of one of these techniques during the. decision step taking place as a 

function of predetermined criteria. 

5. Pecoding 

The invention also relates to the decoding of data encoded according to the 
encoding method described here above. This encoding method is deduced directly 
from the encoding steps- 

Thus, the encoding relies on the reception of a bit stream containing: 
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the de3cnpLion of an initial depiction of the image coming from the 
preliminary encoding (which will be subjected to a symmetrical 
preliminary decoding); 

the values quantified and encoded after DCT transformation 
associated with the selected triangles. 
The weighting coefficients of the matrices may be transmitted in the bit 
stream. However^ preferably, they are known to the decoder. 

The decoding of the values quantified and encoded alter DCT 
ti-ansformation comprises especially the following steps: ; 

Ihe creation of a symmetrical square matrix whose lower part includes 
the dfto.oded cnfifficients of the triangles to be sho\vn, read in the bit 
Stream; 

reverse DCT transformation of the matrix tJius crej^ted; : 
affine transformation of the isosceles rectangular triangle associated 
with the lower part of the matrix towards tlie triangle to bp 
represented. | 
The corresponding decoding especially provides for the reading in the bit 
:>ixe(im received; 

of the number of levels of the hierarchy; 

the identification of the encoding technique used for each of the 
triangles; 

the succession of the differential values of the components associated 
with the nodes of said hierarch]>-a1 meRhing: 

the identification of the arcs on which a diagonal inversion is 
performed. 



19 ; 
APPENDIX 1 

An example of the determining of differential values is shown in figure 4.: 
The figure considers a level 0 triangle whose vertices A, B, C respectively have 
luminance values 210, 150 and 132. This triangle is subdivided into four identical; 
level 1 triangles, thus giving rise to three new vertices E, F and G, located 
respectively at the midpoinLs of the sides AB, AC. and DC of the level 0 triangle; 
The luminance values of the vertices E, F and G. in the image to be encoded ar^ 
respectively 182, 170 and 143. If the luminance values of the vertices E, F, G are 
computed by interpolation of the values of the adjacent vertices, the values 180, 171 
and 141 are obtained. The values introduced into the bit stream are then +2, -1 and 
1 2 corresponding to the difference between the exact values and the interpolated 
values. These values are, however, subjected to preliminarj' quantification to limit 
the number of possible values and thus improve the perfonnance ohaxacteribiics of 
the data compression operation to follow. 

Furthermore, so as not to accumulate quantification errors, the differentia^ 
values are computed from the quantified values of the lower level. 
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APPENDIX 2 

The position is optimized by minimizing a criterion corresponding to the 
luminance divergence between the image to be encoded and the interpolated image.; 
To compute E, we consider points q and 5^ ^ having coordinates x, y and z with zr=Y: 
(luminance component). The criterion E. is defmed as follows: 

where F, is an indicator function eqxial to 1 if q is included between the 
surface formed by the points q of the original image and the surface formed by the? 
points q of the interpolated image, and is zero if not. 

The divergence E. is minimized hy means of aii adaptive pitch gradient 
descent aigoiithm. What has to be done in fact is to find tlie vector X of the points of 
5^ 3 where E is minimal, namely: 

Thb amounts to resolving the following system of non-linear equations: 

This system can be resolved (Meetly by Newton's method. However, it can 
happen that this method does not conveige when the starting point of the iterations is 
far too removed from the solution. It is therefore preferahle to use iterative 
procedures leading to the prodiiction of a series of intermediate positions qi, qs- •■ 
converging towards a local optimum of B. 

This type of iterative procedure runs as follows: the operation to determine 
the optimum position of the peak Sq starts from its initial position <?& , and the 
gradient of E al Ibis point is computed. Since I^^^^ indicates the direction of the 
greatest increase in E, the peak is shifted by a qiaantity otoi^ opposite directioi^, 
and the peak So has the following as a new position: 
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The procedxire thus repeated gives rise to the intarmedialc pusitions qj, q2 ■ • such 
that: 

5 

Advantageously, an adaptive pitch is chosen so as to accelerate the 
convergence. iTie adaptive pitch gradient descent method consists in diminishing (or 
respectively increasing) the pitch when Llic divergence E incrca3e3 (or decreases 

respectively) while complying with a size constraint on the pitch aimed at obtaining 
10 e[an,iijj, cCmax]. The final position q^r is obtained when the diirerence 

between the two successive divergeiit;c values E is smaller then a minimal thresihold 
divergence. The optimum position of the other peaks of the meshing are computed 
similarly. 

As seen here above, the divergence is computed with respect to the real 
15 image. According to an improved version, the divergence E. is computed at each 
meshing level with respect to a specific reference imAge T having similarities in terms 
of frequency content with the meshing level considered. Thus, the basic meshing 
(level 0) represents only the low-j&equency aspect of an imeigc, the reference image 
used to optimize the position of the peaks of this meshing level also has a low- 
20 frequency frequency content. This reference image is obtained by the filtering of the 
real iiiiagc. 

Similarly, with each meshing level, a reference image is associated. The 
frequency content of this reference imcigc is adapted to the meshing considered. The 
reference image associated with the highest meshing level corresponds to the real 
25 image (without filtering). To generate these reference images, an approximation of 
the inf;,nit€ pulse response half-band low-pass filter is used, 

^ ^ — imft : 

L being a sub-sampling factor. 

Thus, a different reference image is used at each meshing level to compute 

30 the value E. ; 

The optimizing of the luminmice values and chrominance values associated 
with the meshmg peaks constitutes another possible improvement of the method of 
the invention. 
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The uplimiziug of (lie valutJis Y, U> V is done through a least error squares 
method aad consists in mininjidng a criterion defined in the field of the image p 
as follows: 



dxdy 



10 



15 



20 



2^5 



Where 

- Sn is a peak with an index n of the multiple meshing, 

- M. is the total ntmiber of peaks of the multiple meshing^ 

I (x, y) represents the luminance value (and respectively the t^hrominance 
value U or V) of the pixel having coordinates (x, y) of the image to be 
encoded^ 

- \1/Sji is the interpolation function associated with the peak Sn? 

- v(Sn) is the optimized Ixmiinance value (and chrominance value 
respectively) associated with the peak Sn- 

If this expression is differentiated^ we obtain the optimised values by 
resolving the following linear system of M equations: 



^(j^y/sr^ix,y), ii/s>.(x,y)^dxdy).v(Sm)=jK2l(x,yy \{/Sm(x,y)Jxdy 



Through the compact base formed by the functions \|fSni, this system of equaiiona can 
also be expressed as follows: 



(1) 



IS 



Where: 

- supp (S^) designates the triangles or mesh units whose vertex or peak; 
the peak Sm? and 

- vei(e) designates the vertices of the triangle e. 
Resolving the above system of equations is equivalent to resolving a matrix system 
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of the foUowmg type: 

AX = B 

Where: 

- A is a symmetric positive-definite matrix^ 

5 - X is a column matrix of the optimized values v(Sxn) with me [L.M]j and 

- B is a coluiim matriTC of the valws of the right-hand term of ihe system (1), 

Since the matrix A is a symmetric positive-definite matrix, it has a unique 
factorization A^LDL"^, L designating a lower triangular matrix with unity diagonal, 
and D being a diagonal matrix for which ail the diagonal coefficients are strictly 
10 positive. The conditioning of the matrix A is furthermore in the neighborhood of the 
unity value. 

The coefficients of L and D can be determined by identifying the coefficients 
located m the lower triangle pan of -foe matrices. 

Besides, we have: 
2Q Lj={Aj-i ^UkWLj^)fDjfor 1<7 < / 

The terms of the matrix X of the optimized values are then determined by the 
formula: 

y-1 M 
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According to an improved embodiment, it can be plaimed to use the profile 
technique to resolve the system AX=B so as to limit the memory size during the 
factorization of A. The profile technique replaces the standard representation of a 
matrix in the form of M X M table by a representation in the form of two vectors. 
Indeed, since the matrix A is hollow (it has a large number of zeros), its 
representation in clasaii: fuixxi is madcquatc. 
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APPENDIX 3 

This diagonal inversion operation pertains only to the convex quadrilaterals. 
Figure 6A. illustrates sltx exemplary non-convex quadrilateral for which the diagonal 
inversion operation is not authorized. 

The diagonal inversion operation on a convex quadrilateral Q formed by two 

triangles Tl and T2 is shown in figure 6B, This operation inverts the diagonal 
representing the edge common to the two triangles Tl and T2. The quadrilateral Q 
then comprises, after inversion, two new triangles T'l and T'2. 

The operation for optimizing the topology then comprises the following 

steps: 

- computing the divergence in luminance E on each convex quadrilateral Q 
formed by two adjacent triangles Tl^ T2 at tlie liighcst mcshmg level, this divergence 
in luminance being equal to the sum of the sums of the luminance divergence of Tl 
and T2, E(Q) - E(T1) + E(T2), 

- inverting the diagonal representing the edge common to the two triangles Ti 
atid T2, so as to form two new triangles T' 1 and T'2, | 

- computing the luminance divergence of the new triangles T^l and T'2 and 
adding these values E(Q) = E(T' 1) + E(T'2), and i 

" keeping^ in the meshing, the two triangles for which the sum of th'C 
luminance divergence values is the lowest. ' 
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CLAIMS 

1 . Image encoding method comprising the following steps: 

- the definition of a hierarchicgi meshing having at least two nested meshing 
levels formed by mesh ixnits defined by penVs (which may be pixels of said image to.. 

5 be encoded) 

- the determining, for each uf said i»esh units, of a piece of error 
information between said image to be encoded and an interpolated image obtained 
from the peaks of the mesh units belonging to the meshing level of the mesh units 
considered; 

30 - putting a stop (33) to the refining of the mesh units having a piece of error 

information below a prcdetcriTdned first threshold; 

- the implementation of a specific encoding (35) for the luesh units having a 
piece of error information above a second predetermined threshold; | 

- the continuation (34) of the refining of the mesh units having a piece o;f 
15 error information above said first predetermined threshold and below said second 

predotcnxdned threshold, ■ 
2. Encoding method according to claim 1, characterised in that said specific 

encoding implements at least one reverj>iblc: txansformatioB> (72, 74). ; 
3* Encoding method according to any of the claims I and 2, characterized in 
20 that said mesh units are triangles (8 1 ). i 
4, Encoding method according to claims 2 and 3, characterized in that said 

specific encoding comprises the following steps; ' 

- the association, witli each of said triangles (SI) to be encoded, of a square 
matrix (84) representing said triangles to be encoded, by means of a first reversible 

25 transformation (72); 

- the application of a second reversible decorrelation transformation (74):to 
each of said square matrices (84) delivering transformed matrices. 

5. Encoding method according to claim 4, characterized in that said step of 
associating a square matrix comprises the foliowmg steps: i 
30 - the affine transformation (82) of a source triangle (81) into an isosceies 

rectangular triangle (83) called a reference triangle; : 

- the creation (86) of a square matrix (84) whose lower part includes dklu 
representing said isosceles rectangular triangle; | 

- the symmetrizing (85) of said square matrix \ 

35 6* Encoding mctliod according to any of the claims 2 to 5, characterized in 
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that said second transformation (74) belongs lo the group compxlslng ; 

- the Karhunen Lo^ve transformation (KLT), • 

- the discrete Fourier transformation (DFT), 

- the discrete (iosine tein^iformation (DCT). 

- and the Walsh-Hadamard transformation (WHT)- 

7. Encodmg method acuuidiug to any of liic claims 1 to 6, characteristed in 
that said specific encoding comprises a step of quantification (75) and of encoding 
the data of the lower part of said transformed matrix. 

8, Encoding method (76) according to any of the claims 1 to 1, wherein at 
least at least one of the following parameters can be modified: 

- said first -flrrcsliold; 

- said second threshol<^ 

- the structure of said nested meshing; 

- a scale factor a enabling an expansion or compression of said square 
matrix in the space domain. 

- the type of qnantification; 

- the quantification pitch. ; 
Encoding method according to any of the claims \ to 8, characterized in 

that said piece of error information corresponds to a difference in luminance, 
representing a squared error or an absolute error between a source triangle and the : 
corrftRpnnding approximate triangle. ; 

10, Encoding method according to any of the claims 1 to 9, characterized in 
that, for each mesh unit^ the following choice is made: : 

- if said piece of error information is below said first threshold, the 
hierarchical division for the me$h unit considered is interrupted (132); \ 

- if said piece of error information is above said first threshold^ but below said 
second threshold^ tiie application of said hierarchical division (134)is continued; 

- if said difference in. luminance \^ ahave said second threshold, said specific 
encoding (135) is implemented, 

11. Encoding method according to any of the claims I to 10, characteri3:ed 'in 
that, for all the meshing levels other than the first level, the data taken into account 
for the encoding are quantified differential values, each corresponding to the 
quantified difference between an optitnized value for the hierarchical level 
considered of the image to be encoded and a value interpolated from the meshing of 
the lower hicrcucliical IcvcL 
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12* Encoding method according to any of the claims 1 to 1 1, comprising a 
step for optimizing the position of the peaks of each of said mesMngs, minimizing a 
predetermined error criterion, 

13, Encoding method according to any of the claims 1 to 12, comprising a 
5 step for defining a quaternary subdivision tree used to identify the mesh units using 

said specific encoding and the mesh units whose contents arc interpolated from the 
values at the nodes defining said mesh units. 

14. Encoding method according to any of the claims 1 to 13, comprising a 
step for matching said mesh units with the contents of the image^ implementing 

10 inversions of diagonals of the quadrilaterals formed by all the pairs of adjacent 
triangles. 

15, A method for the decoding of an image encoded according to an 
encoding method comprising the following steps: 

- the definition of a hierarchical meshing having at least two nested meshing 
15 levels formed by mesh units defined by peaks; 

- the determining, foi* each of sf\id mesh units, of a piece of error 
information between said image to be encoded and an interpolated image obtained 
from the peaks of the mesh units belonging io the meshing level of the mesh units 
considered; 

20 - putting a stop to the refining of the mesh units having a piece of error 

information below a predetennined first threshold; 

- the implementation of a specific encoding for the mesh units having a 
piece of error informalion above a second predetermined threshold; 

- the continuation of the refining of the mesh units having a piece of error 
25 information above said first predetermined threshold and below baid second 

predetermined threshold; 

characterized in that it comprises the following steps: 

- the prclhninary decoding of sdd data encoder! according to said encoding,, 
implementing a hierarchical meshing, enabling the description of an initial 

30 representation; 

- the complementary decoding of said data encoded by means of said specific 
encoding, enablmg tlie refining of said initial representation. 

16. Decoding method according to claim 15, characterized in that said 
complementary decoding step is optional. 

35 17, Decoding method according to any of the claims 15 and 16, charactftrized 
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in that said specific encoding comprisiDLg the following steps: 

- the association, with each of said source triangles to be processed, of a 
square matrix representing said source triangle, by means of a first reversible 
ti"aiisfoniiation; 

- the application of a second reversible decorrelation transformation to each 
of said square matrices delivering transformed matrices, 

characterized in that said complementary decoding advantageously comprises the 
following steps of reconstruction: 

a) the application of a inverse transformation to said second reversible 
transformation on said transformed matrices, delivering said 
reconstructed square matrices; 

b) the association, with each of said reconstructed square matrices, of a 
corresponding reconstructed triangle by means of an inverse affiue 
transformation of said first reversible transfonnation; 

i5 c) the reconstruction of said minimum partition from said reconstructed 

triangle* 

18, Decoding method according to claim 17, characterized in that said square 
matrices are recreated from data of a received bit stream, whu£>e decoded data arc the 
coefficients of the triangle to be reconstructed, which form the lower part of said 

20 matrix. 

19. Decoding method according to any of the claims 15 to 18, characterised 
in that said preliminary decoding provides for the reading, in the bit stream, of at 
leasL uiic of the pieces of information belonging to the group comprisiug: 

the number of levels of the hierarchy; 
25 - the identification of the encoding technique used for each of the triangles, 

the succession of the differential values of the components associated with the 
nodes of said hierarchical meshing 

the identification of die arcs on which a diagonal inversion is made. 
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